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Abstract 

This study uses the current knowledge of marine icing to estimate the amount of spray ice load 
accumulated on the fishing vessel F.V. Destination that sank in the Bering Sea in February 2017. 
The analysis includes wave-impact sea spray and brine-spongy ice models. The input data 
including sea conditions, vessel information and weather information are provided by the US 
NTSB. The main assumptions in the analysis concern the frequency of spray generation, icing 
surfaces and sufficient water delivery for icing. The results of the icing model are strongly 
dependent on the input data provided by the NTSB. Results show that the highest rate of ice 
accumulation is about 2.7 ton/15min which occurred around 1700 AKST on February 10. The 
overall amount of ice accumulated on the vessel, which is calculated using the upper limits of the 
assumptions, is about 154 ton. Using the lower limits of the assumptions, this reduces to 92 ton. 
An analysis was conducted to estimate the effect of the predicted ice accumulation on the vessel 
stability with the result that the upper ranges of the predicted icing create a high probability of 
capsizing due to loss of stability 
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Nomenclature 

c a Specific heat capacity of air (J/kg. K) 

c b Specific heat capacity of brine (J/kg. K) 

c e Equivalent specific heat capacity of brine-spongy ice (J/kg. K) 

q Specific heat capacity of ice (J/kg. K) 

c s Specific heat capacity of brine-spongy ice (J/kg. K) 

e s Standard vapor pressure over the brine surface (Pa) 

H s Significant wave height (m) 

h Convection heat transfer coefficient (W/m 2 . K) 

i Node number (—) 

k a Thermal conductivity of air (W/m. K) 

k b Thermal conductivity of brine (W/m. K) 

k. Thermal conductivity of ice (W/m. K) 

k s Thermal conductivity of brine-spongy ice (W/m. K) 

L Length scale (m) 

L Hb Latent heat of fusion of brine (J/kg) 
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L Hv Latent heat of vaporization of brine (J/kg) 

m Number of the node of solid phase closest to the phase interface (—) 
Nu Nusselt number (—) 

n Number of the discretized grids (—) 

P Wave period (s) 

P a Atmospheric pressure (Pa) 

Pr Prandtl number (—) 

Q e Evaporation heat flux (J/s) 

Q h Convection heat flux (J/s) 

Q r Radiant heat flux (J /s) 

Re Reynolds number (—) 

RH Relative humidity (%) 

S Overall salinity (%o) 

S b Salinity of brine pockets (%o) 

Sc Schmidt number (—) 

T Temperature (°C) 

T a Ambient temperature (°C) 

T b Temperature of brine layer (°C) 

Tbi First interpolated temperature of liquid phase (°C) 

T b2 Second interpolated temperature of liquid phase (°C) 

T f Freezing temperature of brine (°C) 

Tj Initial temperature (°C) 

Tsi First interpolated temperature of solid phase (°C) 

T S 2 Second interpolated temperature of solid phase (°C) 
t Time (s) 

t 0 Initial time (s) 

U r Relative velocity (m/s) 

V Fb Volume fraction of brine (—) 

V a Wind speed (m/s) 

V s Vessel speed (m/s) 

X Phase interface position (m) 

X e The extent of the spray cloud (m) 

a cb Coefficient of the equation of specific heat capacity of brine (—) 

a ci Coefficient of the equation of specific heat capacity of ice (—) 

a kb Coefficient of the equation of thermal conductivity of brine (—) 
a k i Coefficient of the equation of thermal conductivity of ice (—) 
a LH b Coefficient of the equation of latent heat of fusion of brine (—) 
a sb Coefficient of the equation of salinity of brine (—) 

a pb Coefficient of the equation of density of brine (—) 
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a pi Coefficient of the equation of density of ice (—) 

£ Ratio of molecular weight of water and dry air (—) 

£ a Emissivity of air flow (—) 

£ b Emissivity of brine surface (—) 

q a Viscosity of air (kg/m. s) 

p a Density of air (kg/m 3 ) 

p b Density of brine (kg/m 3 ) 

p e Equivalent density of brine-spongy ice (kg/m 3 ) 

pi Density of ice (kg/m 3 ) 

p s Density of brine-spongy ice (kg/m 3 ) 

5 Freezing fraction of solidification (—) 
o Stefan-Boltzmann constant (W/m 2 . K 4 ) 

1. Introduction 

Marine icing poses considerable safety issues such as risks to the stability of marine vessels and 
to the safety of the crew during harsh weather conditions. There are many reports on fishing vessel 
sinkings due to marine icing. Hay (1956) reported the loss of the trawlers “Lorella” and “Roderigo” 
as a result of sea spray icing in the fishing grounds north of Iceland in January 1955. Blackmore 
and Lozowski (1994) reported that this accident was actually the beginning of research on the ship 
icing problems. Shekhtman (1971) reported that 10 Soviet ships in the Bering Sea were lost due to 
instability by the formation of ice in January 1965. Zakrzewski and Lozowski (1988) reported that 
many vessel losses have occurred due to icing. In the 1980s, Jprgensen (1982) reported the risks 
of superstructure icing due to seawater spray and atmospheric icing because of precipitation on 
offshore platforms. Ryerson (2011) covered the hazards of icing, which include: instability and 
loss of integrity of ships and offshore structures, malfunction of the operational equipment and 
communication antennas, slippery handrails, ladders or decks, unusable lifeboats and fire 
equipment, and the blocking of air vents. 

The icing phenomenon on marine vessels in the Arctic and other cold regions: atmospheric 
icing and marine or sea spray icing (Aksyutin, 1979; Shekhtman, 1971). Sea spray icing, or marine 
icing, occurs for two distinct reasons: wave-generated spray and wind-generated spray. Wave 
spray arises from the impingement of waves on marine vessels and structures. Wave-generated 
spray is generally a significant source that is a short and approximately periodic water flux created 
near the bow of vessels and the base of offshore structures. Wind-generated spray is a smaller 
source, but is a constant water flux that is created in the airflow during windy conditions. 
According to previous studies (Zakrzewski, 1986), sea spray is the main reason for icing in cold 
regions. Some parameters have considerable impacts on the icing phenomenon and ice accretion 
during a sea spray event. These include: (1) environmental factors such as wind velocity, air 
temperature, droplets temperature, the freezing temperature of seawater, salinity, the size and 
distribution of droplets, relative humidity, characteristics of swell waves (height, period. 
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propagation direction), sea or ocean conditions, and (2) the characteristics of marine vessels, 
including size and design of the vessels, vessel speed, and the angle between the vessel heading 
and wind/wave direction (Bodaghkhani et al., 2016; Dehghani-Sanij et al., 2017; Dehghani et al., 
2016a, 2016b). 

Ice accretion models for marine icing phenomena require the rate of droplets traveling over a 
marine vessel and hitting freezing surfaces. Velocities and sizes determine the rate of heat transfer, 
vaporization and temperature changes of droplets. Droplet temperatures at the moment of 
impingement on the freezing surface add internal energy to the energy balance for the brine layer. 
In the case of frozen droplets, they add latent heat of fusion to the energy balance equation. The 
rate of mass flux of water to the freezing surface is determined based on droplet sizes and velocities 
(Dehghani et al., 2018; Dehghani et al., 2016b). 

The freezing of a thin layer of salt water and estimating the rate of ice accretion on cold 
substrates are challenging phenomena. Intermittent ice accretion due to wave-impact sea spray is 
a dominant mechanism of ice accretion that depends on heat transfer through the spongy substrate. 
The brine layer starts freezing by releasing and dissipating the latent heat of fusion through the 
substrate. The resultant ice is a two-phase medium including pure ice and trapped brine pockets. 
This type of freezing involves a two-phase substrate called spongy ice and a thin layer of a solution 
of water and salt called the brine layer. Trapped brine pockets, existing between pure ice dendrites, 
fill the porosity of the ice medium and create spongy ice (Dehghani et al., 2018; Schiesser and 
Griffiths, 2009). 

Accumulation of the first layer of spongy ice on a substrate of spongy ice affects the rate of 
ice accretion for the next layers. Estimating the thickness of spongy ice accumulated on marine 
structures has been an important subject of research for many years. The method of the control 
volume, which considers the average of the properties in target systems and is not an accurate 
method for fast phenomena with high gradients of properties, has been the main method of 
approach. A cold substrate of spongy ice is capable of absorbing the released latent heat of fusion 
that can control the rate of ice accretion. A substrate of spongy ice acts as a thermal capacitor in 
the process of ice accretion. In the cooling process, spongy ice becomes colder and the trapped 
brine pockets start to freeze and join the pure ice. The amount of freezing of the brine pockets 
depends on the corresponding cooling process and the final temperature of the spongy ice 
(Dehghani et al., 2017; Dehghani et al., 2018). 

The ice accumulated on the subject fishing vessel is calculated numerically using 
environmental and vessel data provided by NTSB. The numerical scheme uses the current 
published knowledge of marine icing phenomena and a numerical code which is recently 
developed by the authors. The results are strongly dependent on the input data and despite having 
considerable information related to this case, there are some unknown parameters that can cause 
uncertainties. 
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2. Wave-spray ice accretion 

The impact of a wave on a vessel bow creates a thin and unstable sheet of water on the bow surface 
(Dehghani et al., 2016a). The sheet of water splits into several strips. As the strips become unstable, 
they are divided into many droplets. The result is a cloud of spray which is called a wave-impact 
sea spray (Dehghani et al., 2016a). Wave-impact sea spray includes droplets of various sizes and 
velocities. The droplets travel over the marine vessel or the offshore structure to impinge on a 
surface or pass over the platform and return to the sea (Lozowski et al., 2000; Zakrzewski, 1987). 
The spray duration is a short period in the intervals of the wave impact. The wave-impact sea spray 
is an intermittent phenomenon and strongly dependent on time (Kulyakhtin and Tsarau, 2014; 
Ryerson, 1995). Fig. 1 illustrates wave-impact sea spray on a vessel bow. 



Fig. 1. Wave-impact sea spray on a vessel bow (Dehghani et al., 2016a). 

The incoming water on the front edge of marine vessels and offshore structures can be 
estimated based on the Liquid Water Content (LWC) and water flux (Forest et al., 2005; 
Zakrzewski, 1987). The spray cloud starts moving over the platform of vessels or offshore 
structures to reach the cold surfaces. The amount of impinging water droplets on the cold surfaces 
can be predicted by the droplet trajectory method. The duration of the spray cloud is normally 
shorter than the interval between the spray events. This means marine vessels and offshore 
structures are repeatedly subjected to a spray cloud for a short period. Ryerson (1995) reported 
that the average duration of a spray cloud is about 2.7 s. Other studies indicate that the spray 
duration is typically about 2 s (Borisenkov et al., 1975). 

The interval between the spray events can be estimated by assuming spray generation during 
every fourth wave encounter on a vessel bow (Lozowski et al., 2000). In typical sea conditions 
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(Zakrzewski, 1987), the average interval between spray events can be estimated as about 60 s. This 
means that usually in every minute, a spray event of about 2.5 s will occur. For wave-impact sea 
spray, the spray duration is much shorter than the interval periods. The average time ratio of spray 
duration and spray intervals is typically about 1/30 (Zakrzewski, 1987). 

Intermittent spray clouds move over the superstructures of marine vessels and offshore 
structures to reach the cold surfaces and substrates (Dehghani et al., 2016a; Lozowski et ah, 2000; 
Zakrzewski, 1987). The water starts freezing after impinging on the cold surfaces. The freezing 
process is dependent on the atmospheric conditions and the characteristics of the vessels or 
offshore structures. The variables of wave-spray ice accretion are dependent on time. The process 
is almost cyclic and repeats several times. Fig. 2 shows the cycle stages of wave-spray ice 
accretion. This type of ice accretion is different from wind spray ice accretion and atmospheric ice 
accretion which are almost continuous and steady ice accretion processes. 



Incoming spray cloud 
including: 

Spray concentration, 
Droplet itzes and Droplet 
velocities 



Cooling the substrate 
including: 

Conduction heat transfer to the 
substrate. Convection heat 
transfer to the ambient. 
Evaporation heat transfer to the 
ambient and Radiant heat 
transfer to the ambient 



Freezing the brine layer 
including: 

Releasing latent heat of fusion, 
Conduction heat transfer to the 
substrate. Convection heat 
transfer to the ambient. 
Evaporation heat transfer to the 
ambient and Radiant heat 
transfer to the ambient 


Fig. 2. Main stages of a periodic process of ice accretion. 


Incoming spray clouds start the cyclic process by carrying droplets towards cold substrates. 
Droplets impinge on cold substrates and the freezing process starts. A thin layer of brine is formed 
on the cold substrates, freezing starts at the phase interface and the latent heat of fusion is released. 
Heat transfer mechanisms dissipate the released latent heat to the ambient air and substrate. 
Convection, radiation and evaporation heat transfers cause heat flow towards the ambient air and 
heat conduction causes transferring the heat towards the cold substrate. After a short time, the brine 
layer freezes and the new layer of solid ice will be directly in contact with the ambient air. 
Convection, radiation and evaporation heat transfers start cooling down the substrate. The 
substrate will be cold enough to be ready for the start of a new cycle. Wave-spray characteristics 
affect transient ice accretion processes and will be explained in the following sections. 
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The intermittent ice accretion phenomenon can be divided into two individual processes: 
cooling and freezing. Because of the short period of spray events, it can be concluded that for a 
long period, the front surface of the substrate will be dry and the surface will become colder due 
to heat transfer to the environment and sublayers of ice. In this paper, this stage is called the cooling 
stage. In the freezing stage, the outer surface of ice is wet with a thin layer of unfrozen brine. A 
new ice layer, which is confined by the old substrate and the brine layer, exists at this stage. Fig. 
3 illustrates the freezing and cooling stages. 


Wind 

lllllil 

Old Layer of Ice 



Wind 


min 


New Layer oflce 
New Laver of Brine / 






Freezing Period 


Fig. 3. Illustration of the freezing and cooling stages. 


Latent heat released from the phase interface is dissipated during the ice accretion (Shakeel, 
1987; Weeks, 2010; Weeks and Ackley, 1986). During the freezing process, heat conduction, 
convection, radiation, heat transfer due to evaporation and heat released due to solidification are 
relevant heat transfer processes (Kulyakhtin and Tsarau, 2014; Lozowski et al., 2000). The control 
volume around the phase interface includes a thin layer of brine (Horjen, 2015; Kulyakhtin and 
Tsarau, 2014; Lozowski et ah, 2000; Shipilova et ah, 2012). 

In the freezing stage, the heat balance is affected by heat conduction through the substrate of 
ice, heat convection to the ambient air, radiation to the surroundings, heat transfer due to 
evaporation of the brine layer, heat transfer due to the incoming supercool droplets and heat 
released due to the solidification of the brine. It is assumed that the temperature of the brine layer 
is constant and equal to the freezing point (Horjen, 2015; Kulyakhtin and Tsarau, 2014; Lozowski 
et ah, 2000). Previous studies neglected heat conduction through the substrate of ice. They 
considered the other heat flows to be the main contributors in the heat balance (Horjen, 2015; 
Kulyakhtin and Tsarau, 2014; Lozowski et ah, 2000; Shipilova et ah, 2012). Fig. 4 illustrates the 
important factors of the heat transfer in the freezing stage. 
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Fig. 4. Heat transfer factors in the freezing stage. 

In the process of the cooling the outer surface of the substrate of ice, heat transfer factors are 
heat conduction through the substrate, convection and evaporation to the ambient air and radiation 
to the surroundings. Heat flow through the brine-spongy substrate affects the brine pockets and 
the confined pure ice. Cooling will solidify a portion of brine pockets to pure ice and heat flow 
will also melt a portion of pure ice and join it to the brine pockets (Weeks and Ackley, 1986). This 
implies that the cooling and heating cause an exchange of sensible heat and latent heat energy in 
the brine-spongy ice. Fig. 5 shows the heat flows in the cooling stage. 



Fig. 5. Heat transfer factors in the cooling stage. 
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3. Formulation and Modelling 
3.1 Heat transfer 

Models for calculating convective heat transfer for cylinders and planar components were 
presented previously (Horjen, 2015; Kulyakhtin and Tsarau, 2014; Lozowski et al., 2000; 
Shipilova et al., 2012). Important variables in these relations are the relative velocity of air and 
objects, air temperature and the length scales. Evaporative heat transfer, which is caused by 
evaporation of brine, represents another significant variable in the solidification. The factors of 
relative humidity and ambient pressure are needed to calculate the evaporative heat transfer. 
Convection heat flux is calculated using the relevant formula used by previous researchers 
(Kulyakhtin and Tsarau, 2014; Lozowski et al., 2000): 


Q h = h(T b - T a ) 


( 1 ) 


The convection heat transfer coefficient is calculated using: 

h _ Nuk a 
L 


( 2 ) 


Nusselt numbers for cylindrical and planar components are respectively equal to (Szilder and 
Lozowski, 2000): 

Nu = 3.0Re°' 5 (3) 

Nu = 0.036Pr°- 33 Re°- 8 (4) 


Reynolds number of the component is defined as follows: 


Re = HsM (5) 

a a 

Evaporation heat flux is calculated using the following equation (Kulyakhtin and Tsarau, 2014; 
Lozowski et al., 2000): 

Qe = h (if)° 63 1^7(esCTb) - RHe s (T a )) (6) 


Radiant heat flux is calculated using (Lozowski et al., 2000): 


Q r = a(£ b T b 4 -£ a T a 4 ) (7) 

The contribution of brine pockets in the conductive heat transfer through brine-spongy ice 
creates a phenomenon that acts like a thermal capacitor. In a transient condition, heat cannot flow 
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through brine-spongy ice without melting a portion of the ice. Brine-spongy ice cannot become 
colder without freezing a portion of the brine pockets (Shakeel, 1987; Weeks, 2010; Weeks and 
Ackley, 1986). 

Fig. 6 illustrates varying sizes of brine pockets at the end of the freezing and cooling stages. 
At the end of the freezing stage, the brine-spongy ice exists at the highest temperature. The brine 
pockets exist at the maximum volume, and consequently, the volume fraction of brine is at a 
maximum. This stage is a thermally discharged stage. The cooling stage is a stage of charging the 
brine spongy ice thermally. By cooling the brine-spongy ice, a portion of brine pockets starts to 
freeze, and at the end of the cooling stage, the brine pockets exist at the minimum volume. In this 
stage, the volume fraction of brine is minimal. 


End of Freezing Period 


End of Cooling Period 



Fig. 6. Illustration of the role of thermal capacitance in brine-spongy ice accretion. 

Cooled brine-spongy ice has a significant capacity to absorb heat. Due to the contribution of 
latent heat, this capacity is enhanced. The process of cooling brine-spongy ice increases its thermal 
capacity until all of the brine pockets are frozen. The temperature of completely frozen brine 
pockets is about -21 °C (Shakeel, 1987; Weeks, 2010; Weeks and Ackley, 1986). This means brine- 
spongy ice can act as a thermal capacitor in a temperature range from the start of the freezing point 
to about -21 °C. At temperatures less than -21 °C, there are only ice and solid salt phases. At those 
temperatures, the latent heat mechanism cannot be active; therefore, the brine-spongy ice behaves 
like a solid material. 

3.2 Brine-spongy ice 

The governing equation of heat transfer in brine-spongy ice, proposed by Dehghani et al. (2017), 
is obtained as: 


p e c e ^ = V.(k s VT) 


( 8 ) 



( 9 ) 
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This partial differential equation is dependent on the overall salinity and thermal properties of 
brine pockets, ice and brine-spongy ice. Table 1 shows relations which describe the thermal 
properties of ice, brine and brine-spongy ice (Cox and Weeks, 1975; Weeks, 2010; Weeks and 
Ackley, 1986). 

Eq. (8), Eq. (9) and the equations of properties in Table 1 yield 13 equations and 13 unknowns 
leading to a set of linear and nonlinear equations. The unknowns are temperature, brine 
conductivity, brine heat capacity, brine salinity, brine density, brine heat of fusion, ice density, ice 
thermal conductivity, ice heat capacity, spongy ice conductivity, spongy ice density, spongy ice 
heat capacity and volume fraction of brine. Properties of spongy ice are considered as weighted 
properties of ice and brine. The properties in Table 1 are empirical relations found by previous 
investigators (Cox and Weeks, 1975; Weeks, 2010; Weeks and Ackley, 1986). These are 
dependent on the coefficients which are obtained experimentally. Table 2 includes the coefficients 
in the equations of Table 1. With initial conditions, which are the initial distribution of temperature 
in the brine-spongy ice, the boundary conditions, and also the overall salinity of brine-spongy ice, 
this set of 13 equations and 13 unknowns can be solved. The results predict the thermal response 
of brine-spongy ice at each moment of time. 

Table 1. Properties of pure ice, brine, and brine-spongy ice. 


Properties 

Equations 

c b 

« cb o + « cb i T + « cb 2 T 2 + a cb3 T 3 

Ci 

^ciO ^cil T 

C S 

P b X, , P b XT X. 

c b — V Fb + q(l — — V Fb ) 

Ps Ps 

k b 

a kt>o + a k b i T + a kb2 T 2 

kj 

*^ki0 "b *^kil T 

k s 

kbVpb + ki (1-Vpb) 

LhId 

a LH b 0 + a LH b l T + 0CLHt>2 T 2 

s b 

a sbl T + a sb2 T 2 + a sb3 T 3 

P b 

a p b 0 "b ^p b i S b 

Pi 

**pi0 "b ttpii T 

Ps 

P b Vp b + Pi (i — Vp b ) 
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Table 2. Coefficients for the equations of the properties in Table 1. 


Coefficients 



Values 


(0) 

(1) 

(2) 

(3) 

*^cb() 

4211.249 

111.437 

5.125 

93.545 x 10“ 3 

®ci() 

2118.5199 

7.8000 

0 

0 

a kb() 

0.5664 

3.0822X 10 -3 

1.8388X 10 -5 

0 

a ki() 

2.2399 

-10.7517X 10' 

- 3 0 

0 

a LHb() 

-333400.326 

-4958.217 

-29.894 

0 

a sb() 

0 

-17.5730 

-0.381246 

-3.28366 x 10 -3 

a pb() 

1000 

0.8 

0 

0 

a pi() 

917 

-0.1403 

0 

0 


3.3 Wave-impact sea spray 

The movement and frequency of spray clouds determine the rate of water delivery to the cold 
surfaces of the vessel. The extent of the spray is calculated using the formula suggested by 
Lozowski et al. (2000). In this formula, the maximum travel distance of the spray cloud is 
calculated by: 

X e = 2H S + 0.04(1.56P + V s ) 2 — 10 (10) 

Where X e is the extent of the spray cloud, H s is the significant wave height, P is the period of the 
sea wave and V s is the vessel speed. It is assume that one wave among 3 to 5 waves that strike the 
vessel can successfully generate a spray cloud. 

4. Case Definition and Assumptions 

Three categories of information are required to define the case completely for numerical 
simulations. The first category is related to the sea information including significant wave heights, 
sea surface temperature, wave period, wave direction and seawater salinity. The second category 
is the weather condition including wind speed, wind direction, air temperature, relative humidity 
of air and precipitation. The third category of information is the vessel information including vessel 
geometry, vessel speed and vessel direction. The data mentioned above are provided by the US 
National Transportation Safety Board (NTSB). 

The model is defined considering two individual conditions, full load and low load conditions. 
Figs. 7 and 8 show the full and low load conditions. In the full load condition, the vessel is occupied 
by crab cages. For the low load condition, it is assumed that the necessary equipment is on the 
deck and there are no cages on the deck. The geometry of the vessel extracted from Figs. 8 and 9. 
The path of the vessel is shown in Fig. 10. 
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Fig. 6. A picture of the full load condition 



Fig. 7. A picture of the low load condition 



Fig. 8. Side view of the vessel 
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Fig. 10. The path of the vessel 


The significant wave heights and wave period from February 09 to February 11 are shown in 
Fig. 11. The information of the vessel speed and wind speed can be found in Fig. 12. The air 
temperature is plotted in Fig. 13. 
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Fig. 11. The significant wave height and wave period during the vessel’s travel 
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Fig. 12. Vessel speed and wind speed 
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Fig. 13. Air temperature 
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5. Numerical method 

The discretization of the differential equations for the solid and liquid phases for the numerical 
solution is based on a uniformly-spaced fixed-grid method. The freezing fraction, which is the 
ratio of the frozen portion of a grid to the total length of the grid, is defined by 6. The right-hand 
side of the domain represents the brine layer and the left-hand side is the brine-spongy ice. Fig. 14 
illustrates the nodes and grid points from i = 1 to i = n + 1. The phase interface is between i = m 
and i = m + 1. The variable 5, which varies between zero and one, defines the freezing fraction 
in the solidification grid located between i = m and i = m + 1. In order to prevent instabilities 
when the phase interface is very close to the nodes, 5 = 1 or 6 = 0, four interpolated temperatures 
are introduced: T s i> T S 2> Tbi and Tb2- These temperatures provide a smooth variation of the values 
of the discretized equations when the phase interface passes a node. For this purpose, it is assumed 
that temperature varies linearly over the length between two adjacent nodes and the interpolated 
temperatures are calculated based on their distances to the nodes. 



Fig. 14. Discretization of the domain including the brine-spongy substrate and brine layer 


A challenging aspect is discretizing the equations at i = m and i = m + 1, which are affected 
by the moving phase interface. As mentioned, the interpolated temperatures are used as the 
reference temperatures for discretization at i = m and i = m + 1. The freezing fraction and the 
interpolated temperatures are calculated as follows: 

T s i = (1 - SjTm-! + 5T m 
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T s 2 = (1 - 8)T m _ 2 + STjjj.! 

Tbi = (1 — 5)T m+1 + 6T m+2 
T b 2 = (1 — 5)T m+2 + 6T m+3 

dX _ 1 _f k s (T f )+k s (T sl ) T f -T sl _ k b (T f )+k b (T bl ) T bl -T f ~| 

dt Ps(Tf)[l—V Fb (T f )]L Hb (T f ) l 2 Ax 2 Ax i 

It is assumed that the first node, i = 1, which represents the left end of the brine-spongy 
substrate, is kept at a constant temperature. Therefore, the temperature variation at i = 1 will be 
equal to zero: 


(13) 

(14) 

(15) 

(16) 


^i=0 i = 1 (17) 

dt v ' 

Discretization of the partial derivative equations on the interior nodes is obtained as follows: 


dTj _ 1 f k s (Tj +1 ) k s (Tj) Tj +1 T; _^ (TT 

dt p e c e (Tj) 1 Ax Ax ' 


fo( 


ks(Tf)-k s (T sl ) T,-T„ 


- Mt) 


Ax 2 
Tf-2T s i+ T S 2 
Ax 2 


i+Tj-i l 
2 J 


} 


dTj = _ 

dt p e c e (Ti) t Ax Ax 

dTj _ 1 r k b (T bl )-k s (T f ) T bl -T f 

dt Pb( T i) c b( T i) ( Ax Ax 

dTj _1_( k b (Tj +1 )-k b (Tj) Tj F x~Tj , ^-rp Tj F x~2Tj+Tj_x I 

dt Pb( T i) c b( T i) ( Ax Ax '' Ax 2 J 


_ 1_ /'rri N 'I b2 2 I |, j +T| ) 

Ax 2 ] 


i = 2,..., m — 1 

(18) 

i = m 

(19) 

i = m + 1 

(20) 

i = m + 2, ...n 

(21) 


The outer surface of the brine layer is assumed to be an insulated boundary. Therefore, the 
temperature variation will be equal to that of the previous node. 


dTj _ dT n 
dt dt 


i = n + 1 (22) 


The discretized equations, which are arranged based on the MOL, are solved using an iterative 
method by using variable time steps to calculate the temperatures in the brine and the brine-spongy 
ice during the freezing process. 


6. Results of icing calculations 

The mass of ice accumulated on the fishing vessel from 23:15 AKST on Feb. 9, 2017 to 06:15 
AKST on Feb. 11 is shown in Fig. 15. Results vary between Max and Min values corresponding 
to full load and low load conditions. The rate of ice load is shown in Fig. 16. The highest rate of 
ice accumulation is 2.7 ton/15min which is for 1700 AKST on February 10. 
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Fig. 15. Numerical solution of the mass of ice accumulated on the fishing vessel (Destination) 
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Fig. 16. The rate of ice load accumulated on the vessel 
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7. Stability Analysis 

An approximate analysis of the effect of the predicted icing on the vessel stability was conducted 
by estimating the loss of metacentric height (GM) associated with a range of ice accumulation 
values and a range of heights of the effective center of gravity of the accumulated ice. GM is a 
measure of the small angle static stability of a vessel and is a measure of the difference between 
the vertical center of gravity (G) of the vessel and the height of the metacenter (M). Metacenter 
can be thought of as the radius of the shift in center of buoyancy as a vessel inclines and thus the 
distance GM is effectively a measure of the restoring moment arm for small angles of vessel 
inclination. This is the most commonly used index of vessel stability and is relatively easy to 
calculate. 

In the present case, the GM of the vessel is unknown and the analysis calculates a change in 
GM associated with the accumulated ice. The change in GM arises from a shift in G associated 
with the added weight of the ice and from a smaller shift in M associated with the increase in the 
hull submerged volume due to the increased displacement associated with the added weight of ice. 
This approach was taken because the stability and loading condition of the vessel for the departure 
condition are unknown and thus a precise analysis cannot be conducted. In addition certain 
hydrostatic or hull form parameters are unknown (although possibly obtainable) so values typical 
of fishing vessels are used. Obtaining these hull parameters would not materially improve the 
estimate of GM change. 

Fishing vessels typically operate with GM values no lower than 0.35 m and generally not much 
greater than 1,0m. However fishing vessels take on cargo at sea as part of their operation and are 
frequently modified, or loaded without stability assessments. These actions rarely result in 
increasing levels of stability and thus it is not unusual for fishing vessels to operate with GM values 
below the recommended range. In the present case the GM value for the vessel as it existed at the 
time of the sinking or at the time of port departure is unknown. 

Table 3 gives the values used in making the stability loss estimates and the output values for 
a single case. Table 4 gives a summary of output values for all the cases that were considered. 
Figure 17 provides the values of Table 4 in graphical form. 


Table 3. Input and output values for a single Stability Loss Case 


FV Destination - Stability changes based on Icing Predictions 

Assumed vessel condition 

Vessel initial displacement 

650 

tonnes 

Calculated based on parameters below 

634 

m 3 

Calculated volume of displacement 

Assumed initial GM 

1 

m 

Assumed at 1 to calculate deviations 

Estimated KG 

4 

m 

just above main deck (typical value) 

Estimated draft 

3.11 



Estimated KB 

2.07 

m 

2/3 of draft (typical value) 

Estimated Sinkage 

235.2 

t/m 

Calculated based on parameters assumed 
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Inferred BM 

2.93 


= GM-KB+KG 


Assumed Ice accumulation 

Weight of ice 

150 

tonnes 

From Simulation 

Added displacement 

146.3 

m 3 

Calculated based on Sinkage 

VCG of Ice (average) 

9 

m 

WRT Keel 


Changes in vessel condition arising from ice accumulation 

New Draft 

3.75 


Extra weight from ice 

New KG 

4.94 


Change in value from ice accumulation 

New BM 

2.38 


Change in value from ice accumulation 

New KB 

2.39 


Change in value from ice accumulation 

New GM 

-0.17 


Change in value from ice accumulation 

Reduction in GM 

1.17 


Estimated loss of Stability from ice accumulation 


Basic Vessel Parameters 

LOA (overall length) 

33.5 

m 

Taken from Trim and Stability book excerpt 

Beam 

9.8 

m 

Taken from Trim and Stability book excerpt 

Hull Depth (to main deck) 

3.9 

m 

Taken from Trim and Stability book excerpt 

LBP (length for calculation) 

32 

m 

Estimated based on typical values 

Block Coefficient 

0.65 


Estimated based on typical values 

Waterplane coefficient 

0.75 


Estimated based on typical values 


Table 4. Summary of Metacentric height (GM) reductions for all cases considered 


Table of GM Reductions (-m) 

Ice VCG (m) 

9 

7.5 

6 

Ice (tonnes) 

700 tonnes displacement 

150 

1.05 

0.79 

0.52 

125 

0.91 

0.68 

0.46 

100 

0.76 

0.57 

0.38 

75 

0.59 

0.45 

0.30 

Ice (tonnes) 

650 tonnes displacement 

150 

1.17 

0.89 

0.61 

125 

1.01 

0.77 

0.53 

100 

0.84 

0.64 

0.44 

75 

0.66 

0.51 

0.35 


These calculations are estimates based on a range of possible scenarios both in terms of vessel 
condition and in terms of ice accumulation. Some parameters concerning the hull form might be 
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better defined but the loading and operating condition for the vessel at the departure or at the time 
of the loss cannot be known precisely. Since the vessel was known to be heavily loaded at 
departure, two cases of 700 tonnes and 650 tonnes displacement, before icing, were considered. 
Conclusions should be considered as estimates of likely outcomes rather than exact calculations. 

The results of this analysis are expressed in terms of loss of metacentric height for a range of 
vessel conditions and ice accumulations. Predicted losses of GM range from, a low of 0.30 m, for 
75 tonnes of ice accumulation at 6 m VCG, assuming the vessel displacement at 700 tonnes, to a 
high of 1.17 m, for 150 tonnes of ice accumulation at 9 m VCG assuming the vessel displacement 
at 650 tonnes. It appears unlikely that the vessel had an initial GM higher than 0.75m and thus it 
is likely that any ice accumulation greater than 100 tonnes would have placed the vessel in a 
stability compromised condition. Accumulations approaching 150 tonnes would be enough to 
result in a complete loss of stability. 

In the compromised position the vessel would have been vulnerable to capsizing due to 
extreme waves in either beam or following seas and in the case of a complete loss of stability 
capsizing would have been inevitable, regardless of sea state. 


Loss of GM 




0.300 


0.000 


SO IB BO 90 MO 110 130 HO 10D 1V> 160 


Fig. 17. Loss of Metacentric Height (GM) associated with ice load accumulated on the vessel 


8. Conclusion 

The method of active cold surfaces is employed to model ice accretion on FV Destination, a fishing 
vessel that sank in the Bering Sea in February 2017. The data used for this simulation is provided 
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by NTSB and results are strongly dependent on the input data. The main assumptions for this 
simulation are the frequency of the spray generation after wave impact, the area of cold surfaces 
for ice accretion, sufficient sea spray due to harsh weather, transient-rapid freezing of brine and 
interpolation for the gaps between input data points. The simulation shows that the ice accumulated 
on the full load vessel would be between 154 ton and 92 tonnes. Assuming an empty vessel the ice 
load decreases to the range of 77 to 46 ton. The rate of ice accretion before 0600 AKST on February 
10 is negligible. The rate of ice accretion increases to the maximum rate in less than 11 hours. The 
maximum rate of ice accretion occurred at 1700 AKST on February 10, which is about 2.7 
ton/15min. It remains high, more than 1.8 ton /15 min, for the rest of the time. 

Based on the higher range of ice accumulation estimates, the effect of icing on the vessel 
stability, as expressed by the measure of metacentric height (GM), was estimated and it is 
concluded that the predicted ice accumulation would have had sufficient influence to either 
severely compromise, or completely eliminate, the vessel stability. This would result in a high 
probability of capsize associated with the predicted levels of ice accumulation. 
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